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Abstract . In a quasi- two-dimensional p-wave superconductor there find six Cooper 
pairing states which are degenerate within the weak-coupling approach. We show that 
this degeneracy can be lifted by feedback effect favoring the so-called chiral p-wave 
state. This effect is based on the anomalous coupling between charge and current in 
a system with broken time reversal symmetry and parity. 

The discovery of odd-parity Cooper pairing in Sr2Ru04 led to the renewed interest 
in so-called p-wave (spin-triplet) superconductivity. ^ There is a number of other 
systems, such as the heavy Fermion superconductors UPta and UBei3 or the organic 
superconductor (TMTSF)2PF6, where very likely odd-parity pairing is realized. [|[ ^ 
The spin-1 degree of freedom of the spin-triplet Cooper pairs provides a considerably 
wider space of potential pairing states than in the even-parity (spin singlet) case.Q 
The symmetry and effective dimensionality of the electronic band structure plays 
an important role in determining the possible p-wave pairing states. The examples 
listed above cover many of the possibilities: the organic superconductor is quasi-one- 
dimensional, the heavy Fermion compounds are three-dimensional, while Sr2Ru04 
represents the case of quasi-two-dimensional system. 

The large number of possible p-wave states makes their identification for each 
material a difficult task. A simple weak-coupling BCS type of approach can give a 
first guess on the most stable state. Because the condensation energy is in this case 
directly connected with the presence of the energy gap in the quasiparticle spectrum, 
the state with the least nodes in the gap would be most favorable. In one and three 
dimensions the most stable state is unique up to spin rotation. Assuming parabolic 
band structure for the corresponding dimension and using the d-vector notation one 
finds 

d(k) = xfcj, for Id . s 

d(k) = ikx + yky + zkz for 3d 

where the gap matrix is defined as Aj^ = ia20' ■ d(k) and the quasiparticle gap is 
itr[Aj^A]^] = |d(k)p. Obviously the two states are nodeless on the corresponding 
Fermi surfaces. Note that the example for three dimensions corresponds to the Balian- 
Werthammer state or the B-phase of superfluid •^He.Q 

We now consider the case of as quasi-two-dimensional system which is 
characterized by the fact that the Fermi surface is open in one direction, the z-axis. In 
such a system the weak coupling approach does not lead to a unique state, but we can 
find six degenerate states with the same nodeless gap. In Sr2Ru04 their degeneracy 
is lifted by spin-orbit coupling and the corresponding states labeled according to the 
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representation of the tetragonal crystal point group of this compound is given in 
Table | We can distinguish two types of states here: those which have d-vector 
that changes orientation for different points on the Fermi surface belonging to the 
one-dimensional representation Aiu, A2U, Biu and and those which have a fixed 
d-vector orientation but a finite orbital angular momentum, belonging to the two- 
dimensional £'11 representation.!^, ^ Note that the latter is the chiral state, i.e. it 
breaks time reversal symmetry and parity. Since all these states are degenerate in the 
spin rotation symmetric case beyond simple spin rotation transformation, the question 
arises which among them is stable. For Sr2Ru04 experiments suggest the chiral state 
with d II z.|,§ 

Table 1 . Six-fold degenerated states in p-wave pairing symmetry 
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d(k) 


Ai„ 








Bin 


x/cx yky 




Stky -\- ykx 


Eu (chiral states) 


ii{kx i iky ) 



In Sr2Ru04 the loss of spin rotation symmetry by spin-orbit coupling carries the 
main responsibility in picking the stable state. In this letter, however, we assume that 
the spin rotation symmetry is preserved in the normal state so that all states listed in 
Table 1 have the same transition temperature Tc as a solution of the linearized weak 
coupling gap equation. In this case the degeneracy must be lifted in a higher order 
process. A well-known concept introduced by Anderson and Brinkman is the spin 
fluctuation feedback mechanism. |pO|, pl] | If paramagnon exchange plays a dominant 
role in the pairing interaction, the modification of the spin fluctuation spectrum by 
the superconducting condensation also alters the pairing interaction. It was shown 
that this mechanism works in favor of the so-called AMB-state or A-phase in ^He.|^ 
This mechanism applied to the 2D situation turns out to stabilize the time reversal 
breaking state which is indeed the analogue to the A-phase. |jl^, 

Here we would like to introduce an additional feedback mechanism which does not 
exist in a neutral Fermi liquid such as ^He. It is based on the presence of chirality in the 
orbital part of pairing state and we will call it, therefore, chiral feedback mechanism. 
It was shown that in the state d(k) — n{kx ± iky) a Chern-Simons-like term, 

e^j{Aod^Aj + AidjAo) (2) 

appears in the effective low-energy field theory of a static quasi-two-dimensional 
system {i,j — x,y ). ||lj, |l^ Consequently, charge fluctuations generate local magnetic 
field distributions (z-axis oriented) and current fluctuations lead to transverse electric 
fleld distributions, whose orientation depends on the chirality. This property yields 
an additional (anomalous) pairing interaction in the superconducting state which has 
selective power for chirality. This can be seen by the following simple picture. The 
magnetic field induced by the charge of a quasiparticle acts via the Lorentz force on a 
passing- by quasiparticle. ||l^ This force is either attractive or repulsive depending on 
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which side the quasiparticle trajectory is located (Fig. 1). In this way the interaction 

is attractive for one choice of Cooper pair angular momentum (chirality) and repulsive 
for the other. The attractive interaction appears for the same chirality realized in the 
pairs of the condensate. Hence this leads to a positive feedback for the chiral state. 
This effect does not exist for the other states. 




(i) (ii) 

Figure 1 . Two quasiparticles feel (i) attractive or (ii) repulsive interaction depending 
on their relative angular momentum. 



We will now discuss the effect by explicitly calculating its contribution to 

condensation energy immediately below the transition temperature T^- We represent 
the d- vector by d(p) = n^d-yiki/kF where d^i is a complex order parameter and 
kp inverse of the Fermi wavelength. The band structure is simply parabolic ej^ = 
h^{k'^ + ky — kp)/2me without any dispersion along the z-axis leading to a cylindrical 
(open) Fermi surface. We assume that the system has layered structure as Sr2Ru04 
with an interlayer spacing d leading to the density of states, iV(0) = me/27r?i^rf, at 
the Fermi level. The anomalous pairing interaction appearing in the superconducting 
phase is connected with the density-current correlation function which for the 2D 
electrons has the form (in the unit h = c = 1) 

TToj{i'^n,<l)= J (iTd3xe^-"-e^q-^(p(x,T)3, (0,0)) 

x-^ f d^k —(2k +a ) 
= kBT}2 / (0 \3 — ^T>[e^(«t^m+«z^n,k + q)a(iu;™,k)- 

T'^{ii^m + it-n, k + q)T{iiOm-, k)] (j = 1, 2), (3) 
and the Green's functions are 

,^ iwm + e(k) iA(k) 



with £;(p) = ±^e(k)2 + A(k)2 and uJrn = (2m + l)wkBT and i^n = 2mrkBT as the 
fermionic and bosonic Matsubara frequencies, respectively. We express Troj{ii^n,<l) = 
i€ijqif{ivm<^} + VnqjT^{wn,q) due to the gauge invariance. f{iyn,fi) comes from the 
chirality and is written as 

/(^^n,q) = -^e,,^^^^ • (5) 



Close to Tc we can restrict ourselves to the leading contributions in rf^j and obtain, 
, e^fcsTcv^ f d^k -ieijd^id*i,kikv 



where we sum over repeated indices. Note that this expression is only finite for a 
chiral p-wave state. We restrict ourselves to static contributions and comment on the 
dynamical part later. 

We now consider the ^'-dependence. We approximate e(k+q) = e(k) + VF -q where 

vf is the Fermi velocity. In order to evaluate the integral we introduce cylindrical 
coordinates and perform first the integration over the radial part and z-component of 
k, 

., , ie^ksTc , , „ X f d9 kikj 

4me ^ 27r |a;„j|(|a;„|2 + |(vf -qj^j 

where fcj = fci/|k|. For small q we may expand /(q) as 

(8) 



ie2Ar(0)7r 



^C(3) - |c(5)$^«i 



where = <Z^ + <Z^, C = vp/2nkBTc defines the coherence length and ({n) is the 
zeta- function. The behavior of /(q) for S^q± 1 is dominated by the regions of the 
6'-integral for which jvp • q| -C ^nksTc and leads to 

Matching the limiting behaviors together we can approximate /(q) by the following 
form, 

4me {irkBTcY ^\ + 7^2 

which represents the form factor of parity and time reversal symmetry breaking part 
in 7roj(0, q). 

The current-charge density interaction introduced via TToj (0, q) gives an additional 
contribution to the pairing interaction below the superconducting transition. As a 
feedback effect this appears in the GL free energy in a fourth-order correction expressed 
by 

AF{b = kBTcJ ^03i?oo(q)7ro«(q)A, (q)^jo(q), (H) 

following the diagram in Fig. 2. Here Dqq, Dij = 1, 2) is the gauge field propagators 
which in Coulomb gauge are, 



1 d 
^oo(q) = ^— -2 and Ai(q) = (12) 
q + t'pp 



These propagators contain all renormalizations, i.e. Thomas-Fermi screening for the 

scalar potential with the screening length Ztf- Since T « Tc London screening of 
the superconductor can be neglected. We ignore here also the dynamical part for 
simplicity, as it would give the same contributions for all competing states. 




Figure 2. The diagram for the fourth order correction in GL free energy. The 
shadowed circles show the renormahzation by the normal fermionic Green function. 




If we separate the g-integration in q^- and q^-part, we obtain, 



(ff-q^ f dq^ q\ 1 



(13) 



After performing the gz-integration the remaining g_L-integral has a cutoff q±_ ~ l^^. 
This leads to the free energy correction, 



^^8a^re?TF iV(0) . , „ .2 



where we give an expression formahy close to the conventional fourth-order terms in 
order to give a comparison of its magnitude. Here we recover the constants h and c, 
the factor a — /Tic is the fine structure constant and the ratio Tc/Tp indicates the 
strong coupling nature of the correction term, similar to the spin fluctuation feedback 
mechanism. JT]| 

For the chiralp-wave state eijd^id*^ — j2x|A| and zero for all other states (x = ±1 
denotes the chirality). Thus, the correction to the fourth order term is negative definite 
and favors the chiral p-wave state. The ratio between this correction and the usual 
fourth order coefficient is 

Sl3fh o? Tc Itf t.^ 
P ^ nTp d 

which for Sr2Ru04 is of the order 10~^. 

It is easy to see that the dynamical contributions, taking into account Vn 7^ 0, 
does not change the result qualitatively. The corresponding coefficient in the free 
energy, however, increases. We have verified numerically that an increase of one 
order of magnitude is possible. It is clear that other mechanisms, such as the spin 
fluctuation feedback or spin-orbit coupling, would dominate over the chiral feedback 
effect in stabilizing the chiral p-wave state. We would like to emphasize, however, 
that our analysis shows that for a quasi-two-dimensional p-wave superconductor the 
chiral feedback effect, based on the anomalous coupling between charge and current. 
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supports the chiral superconducting phase and, thus, works in the same direction that 
the spin fluctuation feedback mechanism. 
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